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ABSTRACT 


we  show  that  for  multi -dimensional  viscous  flow  computations  the  use  of 
upwind  finite  difference  schemes  can  alter  the  natural  length  scales.  This 
false  scaling  is  related  to,  but  distinct  from,  the  artificial  viscosity 
introduced  by  upwind  schemes.  We  show  that  this  false  scaling  can  account  for 
certain  non-physical  solutions  which  have  been  computed  for  the  driven  cavity 
problem. 
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SIGNIFICANCE  AND  EXPLANATION 


Several  researchers  have  demonstrated  that  upwind  finite  difference 
schemes  for  viscous  flow  computations  give  inaccurate  solutions  when  compared 
with  centered  difference  schemes.  A  major  cause  of  the  inaccuracy  is  the 
artificial  viscosity  introduced  by  the  upwind  differencing.  In  this  paper  I 
show  that  upwind  schemes  can  also  introduce  a  false  scaling  for  multi¬ 
dimensional  problems. 

This  is  therefore  another  reason  to  avoid  upwind  difference  schemes  for 
viscous  flow  computations. 
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UPWIND  DIFFERENCING,  FALSE  SCALING, 

AND  NON-PHYSICAL  SOLUTIONS  TO  THE  DRIVEN  CAVITY  PROBLEM 

John  C.  Strikwerda 

1 .  Introduction 

It  is  widely  known  that  the  use  of  upwind  finite  difference  schemes  for 
equations  describing  viscous  flow  can  introduce  substantial  amounts  of 
artificial  viscosity  at  high  Reynolds  numbers  (see  e.g.  Bozeman  and  Dalton 
(1973),  de  Vahl  Davis  and  Mallinson  (1976)).  The  purpose  of  this  paper  is  to 
show  that  in  multi -dimensional  problems  upwind  differencing  can  also  alter  the 
natural  length  scales  of  the  problem.  In  particular,  in  section  3  we  show  how 
this  false  scaling  can  account  for  certain  non-physical  solutions  which  have 
been  computed  tor  the  driven  cavity  problem. 


2.  False  Scaling 

We  begin  by  considering  a  single  homogeneous  elliptic  equation 


(2.1) 

on  a  rectangular  domain 


32u  32u  3u  3u 


0  <  x  <  1  ,  0  <  y  <  y 

with  u(x,y)  specified  on  the  boundary.  We  assume  that  a  and  b  are 
positive  constants.  The  upwind  difference  scheme  for  (2.1)  is 


ui+ij  -  ^ij  +  ui-ij  t  uij+i  “  *13  *  uij-i 
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which  is  the  sane  as 


(2.3) 
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Now  (2.3)  for  a  fixed  value  of  Ax  can  be  regarded  as  a  central  difference 


approximation  to 
(2.4) 


where 


.2  ♦  82  if”  ,  »  |2  ♦  b  |2 . 0 

8*2  Jy2  *  27 


2  aAx  fl2  bAy 

a  ■  1  +  — r*  and  B  *  1  +  — 


If  we  change  variables  in  (2.4)  by  y*  «  ya/S  we  obtain,  after  dividing  by 


32u  32u  a  3u  b  3u  _ 

<2'5’  Is'TP'sK'w'0  • 

n.  2 

If  we  define  the  Reynolds'  number  of  (2.1)  as  R  •  /a  +  b  then  the 
Reynolds'  number  of  (2.5)  is  R'  ■  /(a/a)2  +  (b/3)2/a2,  and  so  R*  <  R.  This 
is  the  effect  of  artificial  viscosity.  Moreover  the  rectangular  region  for 
(2.1)  has  the  height,  or  aspect  ratio,  of  y  and  that  for  (2.5)  is  ya/0.  We 
describe  this  change  in  aspect  ratio  as  false  scaling. 

Thus  solving  (2.1)  by  upwind  differences  for  given  values  of  Ax  and 
Ay  is  equivalent  to  solving  (2.5)  by  central  differences,  where  (2.5)  has 
both  a  lower  Reynolds'  number  and  different  aspect  ratio  than  (2.1).  Since 
central  differencing  is  second-order  accurate  and  upwind  differencing  is  only 
first-order  accurate  we  claim  that  the  solution  to  (2.2)  for  given  Ax  and 
Ay  is  closer  to  the  solution  of  (2.5)  than  it  is  to  the  solution  of  (2.1). 
This  is  indeed  true  for  the  equivalent  one-dimensinal  problem  for  a  wide  range 


of  parameters  as  is  shown  in  the  appendix. 
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3*  False  Scaling  and  the  Driven  Cavity  Problem. 

We  now  look  at  the  driven  cavity  problem  to  study  the  effect  of  the  false 
scaling  (see  Bozeman  and  Dalton  (1973)  for  a  description  of  the  problem).  The 
equations  in  convective  form  are 


a2»  ,  a2» 

3x2  3y2 


-u 


(3.1) 


32oi  32u  /’dib  3 oi  3i|i  3 o»i 


on  the  square  0  <  x  <  1 ,  0  <  y  <  1 .  R  is  the  Reynolds '  number  and  $  and 
o>  are  the  streamfunction  and  vorticity,  respectively.  The  streamfunction 
and  its  normal  derivative  are  specified  on  the  boundary.  The  top  wall,  at  y 
«  1,  moves  with  unit  speed  to  the  right.  The  other  walls  are  fixed.  Because 
of  the  non-linearity  of  the  sysem  (3.1)  it  is  impossible  to  rigorously  analyze 
the  effect  of  upwind  differencing,  however,  as  a  model  consider 


(3.2) 
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with  a  >  8  >  1.  The  constant  a  is  taken  to  be  greater  than  8  since  it  is 
assumed  that  the  large  velocity  in  the  x-direction  near  the  top  driving  wall 
would  give  a  larger  contribution  to  the  false  diffusion  than  would  the  y- 
components  of  the  velocity. 

As  in  the  previous  section,  let 


and  the  second  equation  in  (3.2)  becomes 


3 


d2w  32cu  R  .9tj>  3u  3ty  3w  . 

which  has  an  effective  Reynolds  number  of  R/afi  and  the  domain  has  an 
effective  aspect  ratio  of  a/0  >  1. 

Now  the  solution  of  equations  (3.1)  for  a  square  driven  cavity  is 
characterized  by  a  single  large  central  vortex  for  any  value  of  R,  Pan  and 
Acrivos  (1967).  Bozeman  and  Dalton  (1973).  de  Vahl  Davis  and  Mallinson  (1976), 
and  Keller  and  Schreiber  (1981).  If  the  aspect  ratio  of  the  cavity  is  greater 
than  about  1.6  the  solution  can  have  (at  least)  two  large  vortices.  Pan  and 
Acrivos  (1967),  Bozeman  and  Dalton  (1973). 

The  use  of  upwind  differencing  for  equations  (3.1)  can,  however,  give 
solutions  which  have  two  large  vortices  for  a  square  cavity  e.g.  Runchel  and 
Wolfshtein  (1969),  Gupta  and  Manohar  (1980),  Shay  (1981),  Bozeman  and  Dalton 
(1973).  This  solution  for  the  square  driven  cavity  is  almost  certainly  r.ot 
correct  as  shown  by  the  careful  studies  of  Bozeman  and  Dalton  (1973),  Keller 
and  Schrieber  (1981),  and  others. 

In  light  of  the  above  analysis,  the  two  vortex  solution  for  the  square 
cavity  can  be  explained  as  the  result  of  false  scaling  which  makes  the 
effective  aspect  ratio  greater  than  1.6.  Indeed,  the  two  vortex  solutions  for 
the  square  cavity  resemble  the  solutions  for  cavities  with  aspect  ratio 
greater  than  1.6  which  have  been  squeezed  onto  a  square. 

It  should  be  pointed  out  tnat  when  upwind  difference  schemes  are  applied 
to  the  divergence  form  of  equation  (3.1)  i.e. 


32u  32u  (3  3^  . 


3  f  3* 


7y  ^ 


the  solutions  exhibit  only  one  large  vortex  for  R  less  than  1000  (Bozeman 
and  Dalton  (1973)).  Why  upwind  differencing  of  the  divergence  form  of  (3.1) 
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should  not  exhibit  the  false  scaling,  but  only  the  false  diffusion,  is  not  at 
all  clear#  It  could  be  that  the  false  diffusion  is  less,  or  that  it  is 
distributed  more  evenly  between  the  two  directions  so  as  not  to  give  a 
noticable  false  scaling# 

4*  Conclusion 

Although  the  analysis  presented  here  is  not  completely  rigorous  it  does 
appear  to  be  useful  in  explaining  the  origin  of  the  particular  non-physical 
solutions  of  the  driven  cavity  problem  that  have  been  obtained  by  upwind 
difference  schemes#  The  analysis  highlights  an  additional  danger  of  using 
upwind  differencing  in  computing  viscous  fluid  flow. 
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Appendix 

Consider  the  one-dimensional  equivalent  of  (2.1), 


(A.  1 ) 


,2 

d  u  du 

— r  +  a  —  -  0 
.2  dx 

dx 


with  0  <  x  <  1,  a  >  0,  and  u(0)  -  1,  u(1)  •  0.  The  upwind  finite 
difference  scheme  for  (A. 1 )  is 


(A. 2) 


v,  -  “i  *  “i-,  .  0it,  -  “i  . 

- -2 -  E - - 


Ax 


with  u0  >  1,  u^i  ■  0  and  Ax  -  1/N.  The  scheme  (A.  2)  is  equivalent  to  the 
central  difference  scheme 


(A.  3) 


<1  +  i~)  I  1  °i“1)  +  a  “i+1  ’  Ui-1 


Ax 


2Ax 


and  (A. 3),  for  fixed  Ax,  can  be  regarded  as  an  approximation  to 
(A.  4) 


2  d2u  du  _ 
°  — 5T  +  a  -r-  -  0 
dx2  ** 


with  o2  m  1  +  — and  u(0)  ■  1,  u(1)  »  0. 


We  will  show  that  for  a  wide  range  of  values  of  a  the  solution  of  (A. 2) 
and  (A. 3)  is  closer  to  the  solution  of  (A. 4)  than  it  is  to  the  solution  of 
(A. 1).  This  serves  to  justify  our  assertions  in  sections  2  and  3. 

The  solution  to  the  difference  equations  (A. 2)  and  (A. 3)  is 
(A. 5)  Uj 

and  the  solutions  to  (A.1)  and  (A. 4)  are 
(A.6) 
and 
(A.  7) 


(1  -  (1  +  aAx)i_N)/(1  -  (1  +  aAx)'N) 


u(x)  -  (1  -  e-a(1-x))/(1  -  e-*) 


u(x)  -  (1  -e-*,(1-x,)/(1  -e"*’) 


respectively,  where  a*  -  a/a  . 
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Zn  Table  Z  we  show  the  values  of  (A. 5),  (A. 6)  and  (A* 7)  for  Ax  «  1/20 
and  several  values  of  the  parameter  a,  at  x  ■  1  -  Ax  -  .95., 

Also  shown  are  the  relative  errors  of  (A.6)  and  (A. 7)  from  (A. 5).  This 
is  an  inverse  error  analysis;  given  the  discrete  solution  (A. 5)  we  wish  to 
know  which  continuous  solution  (A. 6)  or  (A. 7)  is  the  better  continuous 
approximation . 

Mote  that  for  1  <  a  <  100  (A. 7),  the  solution  of  (A. 4),  is  closer  to 

(A. 5)  than  is  (A.6),  the  solution  of  (A. 3)*  On  this  basis  we  justify  our 
claim  of  section  2  that  the  solutions  of  (2.2)  are  closer  to  the  solutions  of 
(2.5)  than  they  are  to  the  solutions  of  (2.1). 

For  a  >  100  the  finite  difference  grid  does  not  have  any  grid  points  in 
the  boundary  layer,  and  it  is  only  due  to  the  simplicity  of  this  example  that 
the  solution  to  the  finite  difference  scheme  is  close  to  the  solution  of  the 
differential  equation.  For  more  difficult  problems  such  grid  spacings  can  not 
be  regarded  as  adequate  since  they  will  not  resolve  any  features  of  the 


boundary  layer 


Table  I 


a 

u19 

From  A. 5 

u( .95) 

From 

rel.  err. 

A.  6 

u(.95) 

From 

rel.  err. 

A.  7 

1.0 

. 07642 

.07715 

1% 

.07642 

0% 

5.0 

.20233 

.22270 

10% 

.20163 

0.3% 

10.0 

.33343 

. 39349 

18% 

.32979 

1.1% 

50.0 

.71429 

.91792 

28% 

.67081 

6.1% 

100. 

.83333 

.99326 

19% 

.76035 

8.8% 

200. 

.90909 

.99995 

10% 

.81112 

11% 
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